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Using  Holder's  niethod  of  frequency  calculation,   the  natural 
frequencies  for  the  first  two  modes  of  torsional  vibration  of  the 
wing  wope  detemined  for  a  representative  conventional  airplane 
(B24-Cy  in  the  oustcmary  manner,  the  fuselage  being  considered  as 
a  rigid  body.     Next,    using  a  method  developed  by  N.   0.   M(yklestad 
of  the  GuggenhciTi  Aeronautics  Laboratory  at  the  California  In- 
stitute of  Technology,   ooaibined  with  ^olzer's  aiethod,   the  natural 
frequencies  for  the  same  two  modes  of  vibration  were  again  deter- 
mined, but  with  the  fuselage  this  tine  being  considered  as  flex- 
ible. 

A  oootnarison  of  results  of  the  two  methods  indicates  that  in 
considering  the  fuselage  as  being  flexible,  a  decrease  in  the 
natural  frequency  of  torsional  vibration  may  be  expected.     For 
the  particular  airplane  selected,  this  decrease  aiiounted  to  6.68^ 
for  the  first  mode  of  vibration  and  to  39.1%  for  the  second. 

Bie  investigation  reported  in  this  paper  was  entirely  theo- 
retical and  was  performed  during  the  1943-1944  school  year  at  the 
Guggenheiai  Aeronautics  Laboratory  at  the  California  Institute  of 
Technology,   Pasadena,  California  under  the  direction  and  super- 
vision of  Dr.   N.   0.  li^klestad,   research  associate  in  aeronautics 
at  the  Institute. 
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INTRODUCTION 

In  the  design  of  aodem  aircraft  for  higher  and  higher  speeds, 
the  designers  are  becoaing  increasingly  aore  interested  in  the 
proble«i8  of  flutter  and  vibration.   One  of  these  problems  is  that 
of  the  torsional  vibration  of  the  wings,  which  is  dependent  upon 
a  number  of  factors,  such  as  {1)   the  taass  distribution  both  span- 
wise  and  chordwisc  of  the  wings  themselves  and  of  all  unite  sup- 
ported either  on  thesn  or  within  them,  (2;  the  torsional  stiffness 
of  the  wings,  (3)  the  torsional  aotaent  anolied  to  the  wings  at  the 
shifting  center  of  nressure  by  the  air  loads,  (4j  the  coupling  be- 
tween the  wings  in  bending  and  the  wings  in  torsion,  («ij  the  tor- 
sioaal  aioment  applied  at  the  root  of  the  wings  by  a  flexing  fuse- 
lage, and  {6)  the  effect  of  conpressibility  as  local  velocities 
over  the  wing  approach  the  velocity  of  sound.   It  is  believed  to 
be  eosEiion  practice  in  the  aircraft  industry  generally  to  consider 
all  but  the  last  two  of  the  factors  enveae rated  above. 

Ihis  laaper  then  has  as  its  objective  the  quantitative  deter- 
mination of  the  effect  on  a  representative  large  airplane  of  the 
fifth  factor  enuraerated  above,  namely,  the  effect  of  the  torsional 
moment  applied  at  the  root  of  the  wings  by  a  flexing  fuselage  on 
the  natural  frequencies  of  the  wing  torsional  vibration. 
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DKFINITTQ-    OF  SY'ABOL^ 

TKin    -  Fuselage  Tiass  concentrated  at  any  station    n,     lbs.    seconds 

squared  per  Inch. 
n      -  Nunber  of  any  wing  or  fuselage  station. 
X     -  fiilass  moment  of  inertia  of  nrlng  in. inch  lbs.    seconds  squared 

or  bending  moeaent  of  inertia  of  fuselage  in  inches  to  the 

fourth  powrer. 
X^    -  C<mvenient  reference  value  of  bending  laament  of  inertis  for 

fuselage  as  a  vrhole. 
E      *  'Modulus  of  elasticity  of  fuselage  bending  material  in  lbs. 

per  square  inch. 
/^„     "•  Angular  deflection  of  wing  at  any  station    r\     in  radians. 
oS      -  Frequency  of  vibration  in  radians   per  second. 
^„      -  f^Miel  length  of  wing  or  fuselage  between  stations    y\      and 

r\  ■\-  \        in  inches. 
5„     -  Shear  at  fuselage  station    n      in  lbs. 
M^    •   Bending  moaient  at  fuselage  station   n      in  inch  lbs. 
Mb    -         "  "         of  fuselage  tail  at  elastic  axis. 

cx^  -  Slope  of  fuselage  axis  at  any  station  n. 

a^,  -       "       "  "  tail  at  elastic  axis. 

oc;  -       "       '•  "  nose     "         "  "     . 

1^^  •  Deflection  of  fuselage  at  any  station  n   in  inches- 

i^b  -      n     M     «     ^g^i   g,^  elastic  axis. 
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(^^  -  Deflection  of  fuselage  nose  at  elastic  axi3. 

y^  -  Change  in  slope  frcfa     r\     to     v^  +  \    due  to  a  unit  force  at     n. 

Vv\  -         "         "       •"  "        n      "       n-f\     "       "  "     "       maaent  at     n. 

dp  -         "         "  deflection  frco     n    to     n+i    due  to  a  unit  force  at   i'^. 

d^^-         "         "  "  ify^ft^  +  ^TTnttw       oioment  at    n 

c^  -  Slooe  of  fuselage  axis  at  extreme  end  of  tail  or  nose. 

^^       i^     -  Coefficients  appearing  in  equation  for       (Xr^  ■ 

M^^   -  Total  coupling  aoment  introduced  into  wing  by  fuselage  at  the 

elastic  axis. 
a„  =  In/^o    "  ^^on-'ii'nensiooal  symbol  for  fuselage  beiiding  maaent  of 

inertia  at  any  station     ri. 

T         — X 
bn  =  <3n+i"''^n"  =-i^^5^== — ==^^        -  Kon-dinensional  aytabol  for  increase  in 

fuselage  bending  moment  of  inertia  fron 

station  n     to  statical     n  +  i  . 

^       -  Variable  distance  froca  station  n    to  any  point  in  panel     X^ 

(between     n     and     n-yi    >. 

K^  -  Torsional  rigidity  of  wing  in  lb.   inches  per  radian. 
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Because  of  the  nature  of  this   investigation  and  the  diffioultica 
involved  in  the  measurement  of  the  effect  of  one  factor  at  a  time 
on  the  torsional  vibration  of  a  wing,  no  experiaental  work  was  uader- 
tfUcen.      Instead,  the  authors  approached  the  problan  frooi  a  purely 
theoretical  viewpoint,   and  the  investigation  was  perfonaed  entirely 
on  that  basis. 

After  a  representative  airplane  for  the  investigation  had  been 
selected,    it  was  necessary  first  to  obtain  the  following  inforaation 
concerning  it : 

A.  The  wing   (Table  Ij,  considering  its  mass  and  the  masses  of 
all  bodies  either  attached  to  it  or  stored  within  it  as 
being  concentrated  at  a  nvraber  of  stations  along  its  span: 

(1)  the  distance  of  each  station  frcm  the  wing  root  in 
inches, 

(2j  the  mass  polar  oionent  of  inertia  X^    about  the  elastic 
axis  of  the  wing  of  the  mass  considered  to  be  concen- 
trated at  each  station  in  lb-inches  seconds   squared, 
and 

(5;  the  rigidity    k^    in  lb-inches  per  radian,   or  its 

reciprocal,   of  the  wing  in  torsion  between  each  eta- 

> 
tion. 

B.  The  fuselage   (Table  II;,   considering  its  mass  and  the  iiasses 
of  all-  bodies  either  attached  to  it  or  stored  within  it  as 
being  concentrated  at  a  nisnber  of  stations  along  its   lengths 


(Ij  the  ddetance  of  each  station  froa  the  fuselage  nose 

in  inches, 
(2;  the  bending  moaent  of  inertia  about  a  horizontal  axis 

perpendicular  to  the  longitudinal  axis  of  the  fuselage 

X^      in  inches  to  the  fourth  pernor,  and 
(S;  the  total  laass   m„    considered  as   concentrated  at  each 

station  in  lbs.    seconds  squared  per  inch. 
After  receipt  of  the   required  inforaation  for  the  wing,  it  was  possible 
to  calculate  the  natural  frequencies  of  the  wing  in  torsion  for  as 
many  modes  of  vibration  as  were  desired,   considering  the  wings  as  being 
built  in  to  a  stiff  fuselage  with  an  extreiaely  high  raanent  of  inertia 
conpared  with  that  of  each  station  along;  the  wing.      This  calculation 
was  actually  carried  out  for  two  iaodes  of  vibration  following  HolEer*s 
-nethod  as   outlined  on  pages  228  and  229  of  T7ef .    1,  an  example  of  which 
has  been  appended  to  this  paper  as   Table  III  with  an  explanation  in- 
cluded in  the  appendix.      The  results  of  this  calculation  have  been 
tabulated  in  Table  'lY  and  plotted  ooa  Fig.   1,  and  show  that  the  natural 
frequencies  for  the  first  two  modes  as  detexnined  by  this  caloulation 
are  53.67  and  71.70  radians   per  second   resnectively. 

T!his  completed  the  first   r:^8e  of  the  investigatior  j  and  with  the 
required  inforaation  for  the  fuselage  then  at  hand,  it  was  possible 
to  proceed  with  Idle  second  phase,  naiaely,  the  oaleulatiem  of  the  natu- 
ral frequencies  of  the  wing  in  torsion  for  8L3  many  modes  of  vibration 
as  were  desired ♦   considering  the  tox^ional  moment  applied  at  the  root 
of  the  wings  by  a  flexing  fuselage.      The  nroblems  immediately  con- 
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fronting!,  the  authors  in  this   phase  of  the  investigation  were  thoae 
of  dete raining   (1>  the  tor3i<»ial  oioaent  oroduced  at  the  root  of  a 
wing  by  a  flexing  fuselage  and   (t)  the  method  of  coupling  this  moaent 
into  the  wing  at  its   root« 

For  the  solution  of  the  first  of  these  probleraa  a  method  devel- 
oped in  "Ref .    2  for  the  ant isjinaie trie  bending  of  wings  was   applied 
to  the  flexing  fuselage,   considering  the  fuselage  to  be  nade  up  of 
two  independent  beaas   extending  in  opposite  directions  fron  the  lo- 
cation of  the  elastic  axis   at  the   root  of  the  wing.      This  niethod  has 
the  advantaee  of  yielding  l^uiediatel^r  the  bending  nfto-aent  at  any  par- 
ticular station  along  a  cantilever  bean  and  the  slope  of  the  beam  at 
that  station  as  linear  functions  of  the  noraal  displaeeaiLent  of  the 
beam.      Consequently,   the  procedure  followed  was,  first,   to  calculate 
the  bending  monents  and  the  elopes,   at  the   location  of  the  elastic 
axis  at  the  root  of  the  wing,  of  both  the  portion  of  the  fuselage 
aft  of  this  location  and  the  portion  of  the  fuselage  forward  of  this 
location  resulting  from  a  unit  downward  displace-ttent  of  the  extreme 
end  of  both  the  tail  and  the  nose.      Ttie  bending  moment  at  the  elastic 
axis  and  the  slope  at  that  location  of  the  after  portion  of  the  fuse- 
lage were  desig-nated  as   Mj,    and    <x^    respectively,  and  of  the  forarard 
portion  of  the  fuselage  as   M^    and   cx.'^,     respectively. 

Hext,   since  the  fuselage  is  actually  a  continuous  structure 
throughout  its  length,  its  slope  on  either  side  of  the  elastic  axis 
must  equal  the  slope  on  the  other  side  of  the  elastic  axis.      This 
leads  to  the  result  that,  since  the  initial  displacements  of  bo^  the 
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tail  end  the  nose  were  taken  to  be  positive  downvrards,   in  order  for 
the  slope  fonrard  of  the  elastic  axis  to  equal  that  aft  of  the  elastic 
axis,  the  slcpe  fonrard  of  the  elastic  axis  «b    taust  be  multiplied  by 
the  ratio        (- ^h/o^b).      This  saae  result  would  have  been  obtained  had 
the  initial  diaplaceraent  of  the  nose  been  multiplied  by  this    ratio 
(■-oCK/oCb)     ;  and  since  the  bendAng  racment  developed  is  a  linear  func- 
tion of  the  displacefuent  of  the  freft  end,  the  bending  nioraent  produced 
at  the  elatstio  axis  b;    the  forward  portion  of  the  fuselage  Mj,   should 
also  be  aulti plied  by  this   sane  ratio,     '".e  then  have  that,  for  the 
continuous  fuselage,  the  bending  moacnts  at  the  elastic  axis  due  to 
the  after  and  forward  portion-^   of  thf^   ^'uselag©  are  given  by  the  ex- 
pre3sions,  Mj,    and  {-<^t>/o(-b)*  ^h         respectively.      However,  these  two 
cocipononts  oppose  one  another?   oonaequently,   in  order  to  detenaine  the 
total  bending  taoaent  M^,   from  the  fuselage  to  be  eoupled  into  the  root 
of  the  wing,   one  must  be  subtracted  from  the  other.      Tf  the  direction 
of  Mb    is  taken  to  be  the  positive  direction,   it  :aay  readily  be  seen 
then  that    M^  =  Mt- (- <^b/«-b)Mb=  Mt,  ^  C»-b/^'v)^\i,.  And  since  the  bend- 

ing jnoncnt  f roai  the  fuselage  at  the  elastic  axis  enters  the  wing  as 
a  torsional  mo:nent,    Mc  i«  the  torsional  aioaent  produced  at  the  root 
of  the  wing,   the  aiiount  of  it  entering  each  side  of  the  wing  being  "2 M^  , 
assuraing  symictrical  twisting  of  the  wing. 

?or  the  method  of  coupling  this  ao-ient   into  the  wing  at  its   root, 
one  side  of  the  wini-  was   considered  as  a  free  body  in  tor?: ion  with  this 
torsional  moraent  of  irM^   aoplied  at  its   root.      Assuming  arbitrary  unit 
angular  deflections  of  the  wing  tip,  the  Holser^s  calculations  ^ade 
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during  the  first  pheise  o£  this   investigation  yielded  th«  follonring  in- 
foxmation  for  each  frequency  selected: 

(1>  The  torsional  moaent  developed  within  the  wing  at  the  first 
station  outboard  of  the  root  due  to  the  rotational  inertia 
forces  within  the  wing  (  ^I„oDA     froa  Table  III>,  and 
(2j  TIae  angle   of  twist  developed  at  the  root  of  the  wing  {/^-, 

froa  fable  III;. 
^ince   in  thia   calculation  the  torsional  moment  developed  at  any 
particular  station  and  the  angle  of  tv/iat  at  that  station  are  given 
as  linear  functions  of  the  arbitraiy  angular  deflection  of  the  wing 
tip,  any  desired  angle  of  twist  at  the  wing  root  can  be  obtained  by 
properly  adjusting  the  arbitrary  angular  deflection  of  the  wing  tin. 
?;ince  the  wint*   can  be  considered  to  be  built-in  to  the  fuselage,   its 
angle  of  twist  at  the   root  should  equal  the  slope  of  the  fuselage  at 
the  wing's  elastic  axis,  and  in  order  to  obtain  this  angle  of  twist 
at  the  root  it  is  neceasary  to  multiply  the  original  arbitrary  anguleur 
deflection  of  the  wing  tip  by  the 'ratio   {^b//^i)'     daving  Jiultiplied 
the  original  arbitrary  angular  "deflection  of  the  wing  tip  by  this 
ratio,    it  is  then  necessary  to  multiply  the  torsional  moment  devel- 
oped T/ithin  the  wing  at  the  first  station  outboard  of  the  wing  root 
by  thia   ratio  also.     Hence,   this  aoraent   is  then  found  to  equal 
(o(.^,//3^XZ^^S}3„     ,   and  adding  this  to  the  torsional  aanent  applied  at 
the   root  of  the  wing  b^'  the  fuselage,   a  re   idual  torsional  aoraent  or 
torquf  on  the -jsrlng  of     M -  aM^  +  (o<b/jfi!,)Xr^«^X  i'  found.      This   residual 

torque  is  then  the  additional  applied  accent  required  in  order  tc  force 
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th©  wing  to  vibrate  at  the  assuned  frequency.      This  value  of  the  re- 
sidual torque  is  then  plotted  against  the  asaunied  frequency)  an"^  the 
process   is  then  repeated  for  other  assisacd  frequencies,    one  point 
on  the  plot  being  obtained  for  each  assuraed  frequency.     A  ocxapletc 
example   of  the  calculation  by  this  niethori   for  one  assuned  frequency 
has  been  appended  to  this   paper  as   Tables  V(a>,  V(b>,  Vl(av,  VI(by 
and  VII, -with  a  brief  explanation  of  them  included  in  the  appendix. 

The   results   of  this   calculation  have  been  tabulated  in  Table  VIII. 

After  a  sufficient  nunber  of  points  have   been  obtained,   a  curve 
may  be  drawn  through  the®  as   has  been  done  in  Fig.    1.     Again  the 
points   at  which  this   curve  crosses  thf?  fTs^qucnoy  axis  detenaine  the 
natural  frequencies  of  torsional  vibration  for  the  wing,   for  at  these 
points  the  residual  torque  becomes   tero,  and  hence  the  additional 
applied  aenent  required  to  force  the  wing  to  vibrate  at  that  frequen- 
cy also  beooaes   sero*     From  lig*    1  it  ean  be  seen  tmat  the  natural 
frequenoies  for  the  first  two  modes  as  determined  by  this   calculation 
are  SI. 42  and  43.65  radians  per  second  respectively,      (^ee  Table  IX 
for  tabulation  of  final   results.;-    These  are  reductions   of  6.68?<d  and 
39.1/0   respectively  fran  the  frequencies   of  the  first  two  raodes   fouhd 
in  the  first  phase  of  this  investigation;     Accordingly,   it  may  be  con- 
cluded that,  whereas  the  consideration  of  a  flexing  fuselage  has  a 
s^all  but  appreciable  effect  on  the  frequency  of  the  first  aode  of 
torsional  vibration  of  the  wing,   it  has  a  very  deciied  effect  in 
lowering  the  frequency  of  the  second  -node.      Probably,   considerinK   the 
trend  of  the  curves  of  Fig.    1,  this   sane  effect  is  carried  on  in  pro- 
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gression  to  subsequent  -nodes  of  vibration;  hence  it  is  the  studied 
opinion  of  the  authors  that  this  effect  should  be  considered  in  the 
calculation  of  the  natural  torsional  frequencies  of  the  wing. 

The  closing  phase  of  this  investigation  was  the  determination  of 
the  fuselage  deflection  curves  lor  each  of  the  two  anodes  of  vibration 
deterained  above.   This  was  accomplished  with  facility  from  the  cal- 
culationa  involved  in  the  detenaination  of  the  bending  motaents  and 
slopes  of  the  forward  and  the  after  portions  of  th«  fuselage  in  the 
second  phase  of  this  investi^^ation.   The  deflection  at  any  station  of 
the  fuselage  is  designated  as  a   ,  and  ooluans  so  headed  may  be 
found  in  both  Tables  V  (a>  and  VI (b>.   Again,  the  values  of  u/ 
given  in  Table  VI  (a>  inust  be  aulti plied  bj   the  ratio  (-cXb/otb)   in 
order  to  ^^i'^e   them  the  correct  magnitude  with  respect  to  those  given 
in  Table  VI (b>.  Fuselage  deflections  are  tabulated  in  Table  X.  A 
plot  of  the  values  of  U-^  calculated  for  the  two  natural  frequencies 
found  in  the  second  phase  of  this  investigation  waa  oiade  and  has  bren 
appended  to  this  report  as  Fig.  2.   A  perusal  of  this  figure  will  in- 
dicate thet  the  deflection  curves  for  the  fuselage  for  the  two  iiodes 

A 

of  wing  torsional  vibration  arf  very  similar,  there  being  no  reflex 
curvatures  along  the  fuaelag?  length  in  either  case.   In  the  first 
mode  the  nose  deflection  is  about  .one-seventh  that  of  the  tail  whereas 
in  the  second  aiode  it  is  almost  twice  that  of  the  tail,  froa  which  it 
can  be  seen  quite  readily  that  for  a  riven  deflection  of  the  nose  the 
fuselage  curvature  will  be  much  greater  for  the  first  mode  than  it 
will  for  the  second. 
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The  relative  angular  deflections  of  the  wing  at  each  station 
along  its  span  can  be  determined  very  readily  by  referring  to  the  col- 
tsms  headed  0     in  Table  III  for  the  calculation  for  a  rigid  fuselage 
and  in  Table  VII  for  the  calculation  for  a  flexing  fuselage.   These 
values  must  be  Jiultiplied  by  the  ratio  {p^\,l(5^    for  each  frequency 
selected  in  the  calculation  for  a  flexing  fuselage,  as  has  already 
been  done  in  the  determination  of  the  residual  torque  acting  on  the 
wing,  in  order  to  detennine  the  actual  nagnitudes  corresponding  to 
a  unit  downward  deflection  of  the  tail.   This  has  been  done  and  the 
results  for  the  two  modes  tabulated  in  Table  XI  and  plotted  in  Fi^.  5. 

Fig,  6  is  a  schematic  illustration  of  the  two  modes  of  vibra- 
tion, assuning  a  unit  downward  deflection  of  the  extreme  tail  in  each 
case. 
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ccncussKms 

In  the  case  of  the  airplane  investigated  herein,  the  consideration 
of  a  flexing  fuselage  has  a  sraall  but  appreciable  effect  on  the  fre- 
quency of  the  first  mode  of  torsional  vibration  of  the  wing,  but  it 
has  a  vejy  decided  effect  in  lowering  the  frequency  of  the  second  mode. 

The  deflection  curres  for  the  fuselage  for  the  first  two  "nodes 
of  wing  torsional  vibratiooa  are  very  siTiilar,  there  being  no  reflex 
curvatures  along,  the  fuselage  length  in  either  case.  However,  for  a 
given  deflection  of  the  nose,  the, fuselage  curvature  will  be  nuch 
greater  for  the  first  mode  than  it  will  be  for  the  seocaid. 

II;  must  be  understood  that  the  a'oove  ocneluaions  apply  only  to 
the  particular  airplane  which  has  been  investigated  herein,   Ihis 
paper  is  not  submitted  with  the  intent  to  show  that  effects  of  ei-iUar 
magnitude  can  be  expected  for  all  airplanes,  but  siaply  that  the  effect 
should  be  investigated  with  the  thought  in  mind  tiiat  it  might  prove 
aporeciable,  oartieularly  in  the  case  of  higher  modes. 
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I.    CiiLCUlATION  FCH   THE  RIGID  TmZlAG?. 

2he  method  used  is  outlined  on  pages  228  and  229  of  !tef .  1,  and 
is  known  as  i=Tolrer*fl  rnethod.  The  wing  data  ("lee  Table  I.;  furnished 
for  the  airplane  in  question  assumed  the  mass  monents  of  inertia  X„ 
of  the  wing  to  be  concentrated  at  7  s^amrise  stations,  the  first  and 
last  stations  being  located  at  the  tip  and  the  root  respectively. 
The  notation  used  is  the  saiie  as  that  for  trie  fuselage  and  is  demon- 
strated in  Fig.   4. 

A  positive   (oliiibing;  pitching  angle  at  the  tip  (    n-i      ;  of 
one   radian  was  assuaed,    (  ^^  =  )     >  and  for  a  given  frequency,  the 
inertia  torque    Xl„tAy3„      was  calculated  for  station    »n  =  \       ,     This 
inertia  torauc  multiplied  by  the  torsitmal  flexibility  for  panel 
length    i,     gave  the  amount  of  twist  or  xhe   reduetioaci  in  angle  /^     be- 
tween stations  n=i       and     n=2     .      This  angle  of  twist  was  then  sub- 
tracted from  /^,     to  give  /^^     ,  the  angular  deflection  at  station 
n=  z       .      Knowing  /^^    ,  the  inertia  torque  at  station    n=2     ,  Z  X^u^/^„ 
was  then  calculated. 

Itie  remainder  of  the  table  was   comnleted  in  like  manner  until 
the   residual  inertia  torque    Zr„(^^„      (at  the  wing  rootj  was  foimd. 
Ihis  value  of  inertia  torque  was  tabulated  in  Table   IV  and  plotted 
against    u)      in  Fig.    1*      The  residual   inertia  torque  is  the  shaking 
moment  which  would  be  requir^^d  at  the  wing  root  to  cause  the  wing  to 
vibrate  toreionall;^'  at  the  assuned  frequency    a)     .  the  case 
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of  a  rigid  fuselage  and  assuaing  syiLaetrio  torsional  vibration,    ^I^J^/3^ 
muot  equal   zero  at  e  natural  torsional  frequency  of  the  wing.      Consequent- 
ly  these  natural  frequencies  can  be  found  by  plotting    ^  X^i^^^   against 
^       to  deteiraine  points  of  intersection  vith  the    u)     axie  as  was  done 
in  Fig.    1.      This  calculation  was  carried  out  for  the  first  two  aiodes, 
the  final  results  aopearing  in  Table  IX  under  "Rigid  Fuselage". 
A  sample  calculation  for     '^TcciS    appears  as  Table  III. 


II.    CALCULATION  FOR  THK  FI£XIBLE  FU5EIAGE 

(a>  Outline  of  Method  Usedt 

The  method  used  here  for  the  fuselage  is  that  derived  in  Pef .   2 
for  the  ant isyiiae trie  bending  of  airplane  wings.     A  brief  resume  of 
the  snethod  follows  herewith. 

Using  the  notation  demonstrated  in  Fig,   4,  the  shear  at  any  sta- 
tion   n      is  given  by 

and  the  bending  moaent  at  any  station    n      is  given  by 
The  slope  at  aixy  station    n+  I        is  given  by 
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and  the  def lecti<ai  at  9iny  station    n  +  \       is  gii'en  by 

where 

Vf„  =    change  in  slope  froa    n+i     to    Ki      due  to  a  unit  force  at    n. 
Vm^=   change  in  slope  froui    n  +  i       to    n      due  to  a  unit  moaont  at  n. 
df^  =   change  in  deflection  from    n-^i      to    n      due  to  a  init  force  at  ri- 
c3m„  =  change  in  deflection  frota   n+i       go    h     due   go  a  unit  iiosient  at  n. 

Tile  'aethod  of  calculation  of  these   oaraaeters  is  outlined  in  sec- 
tion II   (dy  of  this  appendix. 

.  i>ubstiituting  the  expressione  for   S„   and  M„  froa  equations   (1; 
and   (Zj  into  equations   (S>  and  (<«; 

^n.»  -ex,  -  u>\^^Zm,^,   -  J>\m„  Z-^  ^.  ^i  {K-K)  (pJ 

At  the  end  of  the  fuselage,  assume 


80  that      cjt^  =  ^  -  Vp  "I'^m,    -  <^  -  fj.  (Ba> 

Continuing  with  equations    (5 J  and   (6>  in  like  manner  yields 
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wher«  coefficients     "^^  ,    f^      »  <x^     ,  and   k^    are  independent  of    <^      , 
one  eomplete  set  beinc  obtained  for  each  frequency.      The  method  of 
determination  of  these  coefficients   is  outlined  in  section  II   (bj  of 
this  ao-^endix.     For  the  present,   these  coefficients  are  assuaed  to 
be  known  for  any  particular  frequency    u) 

For  antis;naaetric  bending,   the  case  where  the  fuselage  is  being 
shaken  by  a  shaking  aoaent     McostOt    about  the  elastic  axis  of  the 
winr,   the  deflection  at  the  elastic  axis  is   zero    ^^=  o     and  from 
equa-cion   (8^     A  =   1^       .     Tiith.  this  val.ie   of    4*     •   '^^  ©f  the  deflec- 
tions u       aia/  be  found  by  lae&ns  of  equation  (Gy,  as  can  the  bending  ' 
mcmer.t  at  any  station      n     ,    2.  m^oiV  (.x^--)(„)  . 

In  the   ^articular  calculation  with  which  tiiis   paper  is   concerned, 
the  two  quantities  desired  are  the  bending  mooient  caaing  in  frcci  the 
tail  or  nose    M^,   »   and  the  slope  of  the  fuselage  at  the  elastic  axis 

Putting  '11.'*^^^.  "  ^i.  (11^ 

and  m^ui^K    -    K,  (12; 

also  (x.-X,)  -  Z'X,  (13> 


-  21  - 


•b-l  S.W-i  ,  j=b-l  5't>-l 


■men        M,  ='S'[&/Z'ij    -  ''Z'K-^Z  i,l  (U) 

But    't'l&/i'i,l  =  &,a,^l-t-V,)+  S,(l  +  i,-^-i...)  +-&,.,  i,., 

=    Z  [1.    Z  &,]  (16J 

and  «i.nilarly  ZlK<^ZA]  ^  ^'Z  U;  2  K  1 

so  M,  =   *|'  [1.  I  &J    -  <#>*!'  U I  K,l  (16 ; 

deferring  to  Table  VI  (b> 

^  Ks   is  given  by  oolunm  (2>  and 

2.  ^s   is  given  by  col  nan  (6),   the  second  sinaatione  occur- 
ring  In  coluo^is  {Z )   and  {7)  which  colunms  giv* 


I  K 


\,-K  5^^ 


■  K,^Z[i.ZKj  (I7j 

and      Gl=2[i/Z(;j  (18; 


Frcm  this   it  is   seen  that 

Mb  -=     ^'b  -  Kb  "i^  (19  j 

and  froa  oquation  (7;  o^b  -    '^b  ^  -  "^b  (20> 

where    -R^^     and    f^    are  given  on  line   n-=S      under  coluoms    (4j  and   (3 J 
respectively. 

The  .method  was   repeated  for  the  nose  using  the  sane  frequency 
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(Table  Vl!(ay>  to  obtain   M^    aiid     oc^  ,   and  the  four  values  thus   de- 
tei^ined  vere  tabula^bed  at  the  top  of  Table  VII. 

♦Vith  an  assuned  positive   (dosra-wurd;  deflection  of  one  inch  et  the 
tail  the  total  r^oaent  introduced  at  the  elastic  axis  b.    the  fuselage 
is   given  by 

ahi   the  slope  of  th©  fuselage  at  the  elastic  axis  ia   oc^   . 

(b>  T)eterniinatio(n  of  Coefficients 

With  the   original  assuaptions  at  the  end  of  the  fuselage    c<,  =  0 
tLod      u  =  I     ,  frau  equations   (7j  and  (3> 

it  is  obvious  that 

■R,=  I        "^i  =  o  ^.  "  '  k,  =  o        ,   and 

fraQ  equations  (Bay  and  (6a; 

are  obtained  -6^=  i   and  k^=  -^i 

Substituting  equatiaas  (7;  and  (6>  into  equations  (T^;  and  (6>, 

t.,  <?^-^..,=  -^n^  -f.  -  ^Xi"m.(<3.-k.^)  -^N„/Z  r«.(^^-k.<?i)(x.-;(„)     (21; 
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By  eqiaating;  ter^ir,   containing    ^     and  those  not  containing     ^    on  the 
two  aides  of  eRC>i  of  equations    (21;  and   (22 >,  the  following  equations 
arc  obtaineiis 
Tquatiiig  cocjlficieatB  of     ^     : 

Equating  eonstoMt  terniE  x 


Using  substitutions  and  relations  developed  in   (11;,    (12 >,    (13 >,    (15; 
and   {IGji 


X-O-l  i-' 


K,.-K*^f^±K       +    v„„Z[«.ZK.].  (27; 

k.   =  k„    +  iX„   -d.„|K.  -  d^;'£[l.  Z  K.1  (28; 

f„,.  =f„  .v,/Zs^      4    V    2'[l,ZGil  (29; 


n.     -     <^.      -     -i.U      -     d.^    f  ,G.      -    ^M.    Z   U    Z    ^s]  (so; 


All  the  coefficients  -ft^  ,  k.^  »  ^^  »  s^^'^  <^n  ®*^  ^  ;foiaid  by  pro- 
gressive calculation  with  the  aid  of  a  table  (such  as  Tables  VI (a;  and 
VI (b;;  based  on  these  equations. 

(c;  Procedure  Used  in  Filling  cut  Tables  VI (a;  and  VI (b; 


Pfiiranoters    Vp^    ,    v^^   ,     dp^  ,   and    cX^^vtere  calculated   (9««  section 
II (d>  of  this  appendix;  and  written  in  the  spaces   indicated.     Like- 
wise the  values   for   1^    were  entered  in  th<=!   tables. 

Then,  for  a  given  frequency,  values   of  rn„u)     were  calculated  and 
entered  in  the  appropriate  coluran. 

^Teart,  in  line    r\=\        ,  the  following  values  were  entered  in 
colurnns    (1^    (4;  and   (*5>  respectively i 

k.  -=  o  ^^=  \  ^i"  ^ 

and  in  line    n^^z     under  oolunm   (l^j  k^=  ^, 

Ihc  table  was  then  worked  aoroGS  frora  left  to  right  starting  with 
line    r\-\      then  proceeding  with  line    n^z      etc.,   each  step  being  in- 
dicated in  the  colunn  heading. 

In  computing  values  to  enter  in  col'jnms    (1>,    (5/,    (5j  and   (7 ),   on« 
must  rcmcrabcr  to  use  information  appearing,  in  the  preceeding  line. 

The   remainder  of  the  steps  are   self  exolanatory ,   the  desired 
quantities  of  the  calculation  bein;^  M^,    and    <x^  . 

(dj  Calculation  of  Parametem   Vp      ,    v,^     ,  dp      ,  and    6^   . 

The  fuselage  data  received   (Table  II.;  indicated  bending  mcpaents 
of  inertia  equal  to  eero  at  each  end  of   the  fuselage,   but   in  order 
to  more  nearly  approxltiate  the  probable  raofnent  of  inertia  distribution 
in  the  regions  fron  the  extreme  ends   to  the  next  stations  inbcard, 
a  trarjesoidal  distribution  over  these  regions  was  assumed  in  both 
cases  with  the  end  ordinates  approximately  half  the  value  ol  the  next 
ordinates  inboard.      Ihe  assuaed  values  were 

•I.     (TAIL)     =    50Q>    In/^  I.  (NOSE")  =    3,500    in'*- 
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A  fuselage  stiffness  cunre   (suoh  as  Fig.    5>  could  be  drawn  for 
the  given  and  assumed   (endy  value e  of  binding  motixent  of  inertia,  where 
the  bending  jooment  of  inertia  is  plotter!!   against  fuselage  distance 
as  abscissa,   such  that  the  areas  between  aucceeding  stations  would  be 
trapezoids.      If    I^    is  taken  as  a  convenient  refe^^nce  value  of  the 
bending  umient  of  inertia  for  the  fuselage  as  a  whole,   then  at  any 
point  between  stations    n       and     n+i 

T'sin^^-  the  mocnent  area  laethod: 

0 

Heferrinf  to  Tables  V(a>  and  (b>,  EX^  was  taken  as  to'**  •  Coliriins 
(1)  and  {2 J  were  filled  in  with  values  of  1„  end  o^  fraa  the  data 
furnished.      The  value   (a^fti,)    in  column   ($ )  is  determined  from 

Due  to  lack  of  availability  of  a  seven  olacc  table  of  natural 
logarithns,-  cowaon  logarithms  were   used  and  values   converted  to 
natural  logarithms  in  coliams   (8>  and   (9j  by  the  relation 

The   remainder  of  the  table  is   self  explanatory  and  follows   frocn 
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equationa   (Sl>,    (52;  and   (35;.      The  desired  parscaeters  appear  in  col- 
mna  (G;,    {12J  and  (18 j. 

(©;  'Procedure  ^Taed  in  Filling  out  Table  VII. 

/^ain  ur>lng  the  naiie  frequency  as  was   used   in  Tables  VI (a;  and 
VI (b;,  HolEor's   calculation   (as  explained  in  section  I  of  this  appen-  ■ 
dix;  was   repeated  in  Table  VII,   the  desired  quantities  being  the 
angular  deflection  of  the  T^ing  at  the   root  /3-,    and  the  inertia  torque 
at  the  next  station  outboard  frcn  the   root     '^X^'^/S^  Tbiis  inertia 

torque  adjustc^d   so  as  to  i;ia>e  the   anarular  deflection  of  thr  ^ring  at 
the   root  equal  to  the   slope  of  the  fuselare  at  the  claistic  axis   is 

(Zl„ti>X}  (**^''/^-')  *  ®^^  this   adjusted   inr-tia  torque  aided  to  half 

the  bend^n;;  -nn^^^-nt  introduced  at  the  elastic  axis   oy  the  fuselage 

i"  ^c      gives   the  residual   torque  actinp   at  the  'winn;   root  which  would 

be   required  to  ^ako  the  winr  vibrate  torr,ionally  at  the  chosen  frequency. 

At  a  natural  frequency  of  the  system  this   residual  torque  is   zero. 
Table  VIII  gives  a  tabulation  of  the  results  of  this   calculation. 
Vig.    1   shows  a  T>lot  of  this    reaidt^il   torque  against  c^       and  the  natural 
frequencies    (first  and   second  rnodes;  occur  when  this   curve   crosses  the 
60        axis.     The  final  results  are  tabulated  in  Table  LX  under  "Flexible 
Fuselage". 
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MODE 

.RIGID 
FUSELAGE 

Ci) 
fRADIANS/SEC) 

FLEXIBLE 
FUSELAGE 

0) 
(rad<ans/sec) 

PERCENTAGE 
DECREASE 

FIRST 

33.472  15 

32.42148 

G.G8 

SECOND 

71.7  0349 

43.(&5\7  1 

39.  1 
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NUMBER 
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u)=43.G»4<i3057 
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FOR.    TAIU 
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06,84464- 

I.  240J»-38 

3 

151.00 

0.4  4^298 

0.06465^ 

0.3S,3O22     . 

0693995 

4 

225.00 

0.206  II  Z 

0.03042) 

0.154  186 

0.2793(39 

5 

' 

258.56 

0.117336 

+  O.OI7i5Z 

0.0787e» 

+  0.14274-3 

Q 

9 

305.00 

0 

0 

0 

0 

8 

349.50 

-  a  OOO047 

-  0.062598 
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5 

554.00 

0.377  75-3 
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0.G2C.639 

0.SO735'<6 
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0.7-/70  1  1 

0.63&222 
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75&.O0 
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0.8373<i»2 
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791.50 

l.OOOOOO 

l.OOOOOO 

TABLL  X  -  FUSELAGE  DEFLFCTIOKS 


-   40  - 


NUMBER 
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FIG.    4  -  SKETCHES  SHOEING  NOTATION  ITJED   IN  BENDIMG  CALCULATIONS 
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FIG.    5  -   SA-iPLE  FUSELAGE  STIFFNESS   CUP^^'E  AND  NOTATION  USED 


-  45  - 


o 


o 

-a 

M 


(r* 


i9 


^ 


VMS  DOB 

' 

-- 

« 

-«•- 

■ 

UN       2 

596 


Tue 


_  J   m'''  Furer  15483 

The  quantitative  effect  of  a 

f   JTJorsional  modes  of  the  wings  of  a      '^ 
^fllar^  airplane        _  ss  oi   a 


